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Abstract .Ba,sing on the thickness (profile) function, previously obtained for the realistic Fermi type distribution 
of nucleons in nuclei, calculations are made of the microscopic eikonal phases of the nucleus-nucleus scattering 
and the total reaction cross- sections. In so doing, the phase is deduced to the one-dimensional integral provided 
CM that the Gaussian density distribution for the projectile nucleus and an arbitrary shape of the thickness density 

for the target nucleus are used. The problems of obtaining parameters of the "point" nucleon density are 
considered. A possibility of approximating the realistic densities by the " surface- matched" Gaussian functions 
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and the dependence of cross-sections on the nucleon-nuclcon interaction radius are discussed. The in-medium 
effects and the role of the trajectory distortion are studied. Conclusions are made on physics of the process, 
and comparison with experimental data is made with cross sections calculated by using the developed method 
where no free parameters are introduced. 
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oo '■ 1 Introduction 
o 

In nuclear physics, the Glauber-Sitenko approach [1,2] is used with regard to several its modi- 
fications for investigation of nucleus-nucleus collisions at energies of the order of 10-rl00 MeV 
per nucleon of an incident nucleus. In this case, one can determine the eikonal phase by both 
introducing a model optical potential and using the microscopic approach where it is expressed 
through density distributions of nuclei and the nucleon-nucleon scattering amplitude. 

The parameters of the phenomenological potential U opt (r) = V(r) +iW(r) are usually fitted 
by comparison of calculations with the experimental differential elastic scattering and total 
reaction cross-sections. Nevertheless, the problem of ambiguities of obtained parameters still 
remains [3]. For example, it was shown in jlj that the total reaction cross-section obtained 
analytically for the Woods-Saxon potential depends, in principle, on two combinations RWq 
and R/a of three parameters, the radius R, diffuseness a and the power Wo. This fact enables 
one to specify these latter quantities in the wide limits of their available values. As to the 
microscopic approach, it does not include, in principle, free parameters, and enables one to 
calculate the eikonal phases themselves rather than potentials of scattering. It allows one to 
make predictions of the total reaction cross-sections, including, e.g., those with participation of 
radioactive nuclei, the important problem related to transmutation of radioactive waste E] . 

Beginning with the earlier work [7], the microscopic approach was applied for calculations 
of proton-nucleus cross-sections. Later on, based on the respective theory of the multiple 
scattering of nucleons by nuclei [TJ |2] , it was generalized in [HI E| to the processes of nucleus- 
nucleus scattering. This approach was often employed for analysis of interactions of light nuclei 
with nuclear targets which led, in particular, to the discovery of the neutron hallo in 6 He, 9 Li, 
the proton hallo in 11 Be, and the establishment of nuclei with excess of neutrons or protons 
(see, for instance, [TU] . and the review papers [TTJ H2|). Much attention was also given to 
studying the mechanism of scattering of nuclei and, in particular, to effects of deflection on the 
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straight-line trajectory of motion E|. The role of higher order corrections to the eikonal 
phases [HI] and effects of nuclear clusterization ^1 were examined, too. 

In the majority of such works the use of Gaussian functions (G-functions) for nuclear den- 
sities is typical because they make it possible to separate variables in the multidimensional 
integrals for eikonal phases and, therefore, to obtain results in the analytic form. This is 
the main reason why G-functions and their modifications are used in calculations of cross- 
sections even for heavy ions although in this case, of course, the physical reasons require 
taking the functions of extended shape. In nuclear physics, the Fermi function (F-function) 
Uf(r) = [1 + exp(r — R/a)]~ l is usually used as the most realistic function for description of 
densities and potentials of medium and heavy nuclei. Unfortunately, in this case the analytic 
calculations encounter difficulties. For example, it is impossible to separate variables in the 
same multidimensional integrals for phases. Nevertheless, F-functions are applied not only in 
numerical, but also in analytical calculations. In the case of heavy ions these functions are es- 
pecially needed since they are most realistic ones for description of the shape of potentials and 
densities in the periphery of collisions of nuclei, the region that forms the elastic differential and 
total reaction cross-sections. With reference to microscopic approaches the problem is posed to 
develop analytical methods for calculating cross-sections when using the realistic densities for 
nuclei with their individual parameters being known from other experiments. In so doing, one 
can predict with confidence the cross-sections for different combinations of colliding nuclei, and 
thus to study the genuine mechanism of their interaction. Just such is the goal of the present 
work. 

In Chapter 2 presented is a series of modifications of the origin (initial) microscopic formula 
for a scattering phase. This is important for both understanding the mechanism of nucleus- 
nucleus scattering and calculations free of assumptions carried out in many works. In Chapter 
3, the explicit expressions for the so-called profile functions of nuclear densities are given in the 
form of the Gaussian, uniform, and symmetrized (SF) Fermi functions. It is shown how one 
can reduce the 4-dimensional convolution integral to the one- dimensional one if the density of 
incident (light) nuclei is presented in the form of G-function and the density of a target nucleus 
in the form of SF-function. For arbitrary forms of densities the phase is obtained as a one- 
dimensional integral with the Fourier-Bessel transforms of profiles of densities. The explicit 
form of such a transform is given for the profile of the SF-density. An inference about the 
applicability of the so-called modified Gaussian densities is drawn. Chapter 4, is devoted to 
the problems of using nuclear densities obtained from nuclear form-factors in electron-nuclear 
scattering, the effects on cross-sections of a choice of the NN-force radius, the distortion of the 
trajectory of scattering, and the in-medium effects. In Chapter 5, we give the comparison with 
experimental data and general conclusions. 

2 Basic formulae 

In the framework of the eikonal approximation and the microscopic approach one can obtain 
the total reaction cross-sections as follows [HI2] : 
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where the phase 



X (b) = a NN l{b) (2.2) 
is determined by the isospin average total nucleon-nucleon cross-section 

_ Z p Z t a pp + N p N t a nn + {Z p N t + N p Z t ) 

^np ( c\ q\ 

aNN ~ A^A t (2 ' 3) 

and the folding integral that, in the case of the nucleus-nucleus scattering, has the form [HI HI] 

1(b) = J d 2 s p d 2 s t p° p (s p ) p° t (s t ) /(£), e = b - s p + s t . (2.4) 

Here vectors s, £ are displayed in the impact parameter b plane perpendicular to the oz axis 
along the momentum kj of the incident nucleus 1 , and p°(s) are the so-called thickness (profile) 
functions of density distributions of centers of nucleons ("point-like nucleons") of the incident 
nucleus (with the atomic number A p ) and the target nucleus A t 2 . The thickness densities are 



given as 



p(s) = dz p(Vs 2 + z 2 ). (2.5) 



The point nucleon profiles p°(r) differ from the matter distributions p{r) in nuclei composed of 
real, i.e. "dressed" nucleons. So, when using the convolution integral ()2.4j) we need be concerned 
with obtaining the point-like densities p°(r) from the "experimental" nuclear densities p(r). 
Just for p(r) a large set of tabulated data exists obtained from electron-nucleus scattering 
data 3 , and our goal is to develop the approach in such a way that to use in calculations the 
table data for p(r) and do not introduce free parameters. 

The function /(£) determines the form of the nucleon-nucleon interaction amplitude 

/(0 = (V^a N y 2 e-£V4, 4 = \r 2 Nrms . (2.6) 

Here r 2 Nrms is the root-mean-square radius of NN-interaction, and a 2 N = 2(3 is expressed by 
the shape parameter (3 of the scattering amplitude of nucleons 4 in the form exp(— q 2 f3/2). The 
values of (3 at energy of the order of 1 GeV are in the interval 0.21-H0.23 fm -2 JH] which means 
that r 2 Nrrns = 0.63 -hO. 69 fm 2 . In our case the nucleon-nucleon forces act in nuclear medium, and 



1 I n |17| a similar expression was obtained in the model of interacting tubes of the nucleon fluxes in collided 
nuclei. 

2 Here s and r are vectors in the 2- and 3-dimensional spaces, and r 2 = s 2 + z 2 . Then, p(r) and p°(r) are 
the density distributions, and p(s) and p°(s) are their profiles, respectively. 

3 In general, one-particle densities p(r) depend on coordinates in the respective center-of-mass frame of a 
nucleus. However, as usual, in analysis of form factors of nuclei, one omits the respective factor of the center- 
of-mass motion exp(q 2 (q 2 ) /6A), where (r 2 ) and A are the root-mean-square radius and a mass number of the 
nucleus. Therefore, the tabulated p(r) appear to be distributions of a nuclear charge (or matter) in the field of 
a fixed nuclear potential. At small q and large A densities in both systems coincide. 

4 The amplitude is /jv(?) = /iv(0)/(?), where /jv(0) = (fc/v /47r)CTjvjv(* + &nn), and fcjv is the rel- 
ative momentum of colliding nucleons. For f(q) = exp(— q 2 a 2 N /A) the Bessel- Fourier transform /(£) = 
(2tt)~ 2 J exp(—iq£)f(q)d 2 q follows ea. (|2.6|l . and in the zero-range approximation (ajv = 0), when f(q) = 1, 
one gets the delta-function in 2-dimensional space, i.e. /(£) = <5( 2 )(£). 
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to take into account their influence the correction factor f m is additionally introduced under 
the integral. Later on we will touch this problem in detail. 

The convolution integral (|2.4J1 has a similar form as the 6-dimensional "double folding" 
integral in calculation of the nucleus-nucleus potential in both the cases we are led to 

search for the ways of separating variables in integrands. In Chapter 3 we show that the integral 
(|2.4j) can be calculated explicitly if both densities are the G-functions, or it is reduced to the 
one-dimensional integral when one of the densities has the Gaussian form. At the same time, 
there exists a standard procedure to transform such integrals to one- dimensional ones in the 
momentum space. To this end, in each function in (|2.4|) one should perform the two-dimensional 
Fourier-Bessel transformation 



u (s) = —}— f e ~ iks u(k)d 2 k = / J (ks)u(k)kdk, (2.7) 



where 



Then f|2.4|) becomes 



where 



(2vr) 2 J v ' 2tt 



u{k) = I e iks u(s)d 2 s = 2n I J (ks) u(s)sds. (2.8) 
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Tib) = — I kdk J (kb) p%(k) p%(k) f(k), (2.9) 



f(k) = exp (-k 2 r 2 Nrms /6) . (2.10) 
Next, using the convolution formula for the nuclear thickness density 

Pi(s) = J d 2 s N p N {s N ) p°(\s - sjv|), (2.11) 

where pn(sn) is the nucleon thickness density, one obtains with the help of ()2.7|) the following 
result: 

p i (k)=p N (k)p\(k). (2.12) 
For the Gaussian density of a nucleon with the rms radius one has 

p N {k) = exp(-^|^i). (2.13) 

Then, ()2.9j) results in 

oo 

1(b) = —J kdkJ (kb)p° p (k)p t (k) exp( (2.14) 

o 

-2 „2 2 



T 



r Nrms r 0rms- (2-15) 



If it is considered that r 2 Nrms and rQ rms coincide, then r 2 = and so 

oo 

1 



1(b) = — kdk J (kb) p° p (k) p t (k), (2.16) 
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and in the coordinate representation 

oo 

1(b) = J d 2 s P ;(\b-s\)p t (s). (2.17) 



As a result, we have obtained the expressions for convolution integrals (j2.14j) . (j2.16j) . (j2.17j) . 
where the thickness functions pt{s) (or pt(k)) come into being instead of profiles of the target- 
nucleus "point-like" densities. Just for them the respective "experimental" densities are usually 
known from tables where they are parametrized for heavy and middle-weight nuclei in the form 
of a Fermi function. In principle, one can employ eqs. (j2.12|) and (j2.13|) for incident nuclei 
as well, i.e. one uses p° p = p p /pN- But then under the integrals (|2.14j) . (j2.16j) an increasing 
Gaussian function appears, and if for the profiles of both the densities the functions with 
realistic exponential asymptotics are taken, then the integrals will diverge at the upper limit. 
Of course, one can act formally, i.e. either "cut" integration at the point where an integrand 
starts to increase, or change the Gauss-like nucleon form-factor (j2.13|) by a dipole formula (see 
below eq. (|4.2j0 . On the other hand, if for the density of one of the nuclei one takes a Gauss 
function, then no divergence arises. 

Finally, it should be mentioned that in some papers the so-called zero-range approximation 
( r Nrms = 0) i s used. This leads to the convolution integral in the form (j2.9J) with f(k)—l, 
or, in the spatial coordinates, in the form f)2. 17j) with pi (s) in (j2.17|) stead of pt(s). A rougher 
approach is when both the densities in (J2.17|) are supposed to be the nuclear ones p(s). Now 
we see that such approaches are not necessary to be used, and they themselves distort the true 
mechanism of nucleus-nucleus scattering. 



3 Eikonal phases for realistic density distributions 

As mentioned above, to obtain analytic expressions for phases and cross-sections the Gaussian 
(G-functions) density distributions and their profiles are used 

Pc(r) = p G (0) e- r V4, PG (0) = A/(^a G )\ (3.1) 

p G (s) = (V^a G ) p G (0) e- s V4, 4 = l R l mgi (3.2) 

where the only parameter a G is determined by the root-mean-square nuclear radius -R rms 5 . 
We also list the functions of uniform density distribution and the relevant profiles 

Pu(r) = p u (0) B(R U - r), p u (0) = 3A/4rrRl (3.3) 

p u (s) = Pu (0)^Rj^ Q(R U - s ), Rl = ^R 2 rms , (3.4) 

which are sometimes used for middle and heavy nuclei. 

In principle, the realistic density having the Fermi distribution (F-function) can be approx- 
imated by a sum of Gaussian functions with by fitting expansion coefficients and parameters 

5 For the point density p° G (r), the parameter a^. 2 = \V? rms can be expressed through the respective rms- 

rnrlinQ 7? 2 - f? 2 t- 2 

lauiun i\* rms — K-rms 'Orras 1 
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a G . Such a procedure was suggested in [213, an d i n EI! this one was proposed not to densities 
but directly to profiles of Fermi functions. Unfortunately, this procedure must be repeated 
at every new set of parameters R and a of the initial F-functions. However, one can remind 
that for heavy ions, both differential elastic scattering and total reaction cross-sections are first 
determined by the behavior of phases in a periphery of collisions. It seems likely that the first 
attempt to model the tail of the F-distribution Pf{t) with the help of one Gaussian function 
was made in [22] ■ Later on in |23| . not a Fermi density but its profile Pf(s) was reproduced in 
such a way. For this aim the G- and F-functions were matched in their periphery to find two 
parameters a G p G (0) of the so-called modified G-function 

Pais) = (V^clq ) p e (0) e- fi2 /4. (3.5) 

This function is not normalized, since its parameters are no longer connected mutually in such 
a way as in (J3.1|) . In the general case, when matching po{s) with an arbitrary form of the 
extended function p(s) at two points S\ and S2 we obtain two parameters of (J3.5j) 

Pg(0) = (V^a Q y p(si) exp(sl/a%), (3.6) 



a G = ] — r~\ — / ( \ ■ ( 3 - 7 ) 

Llnp(si) - lnp(s 2 )J 

In j2H] , the points of matching the G- and F-profiles have been taken as s\ = c and S2 = c + Ad 
where c and d are the radius and diffuseness parameters, and the profiles p° F (s) were obtained 
by numerical integration in (J2.5|) . Herewith d was assumed to be the same d = 0.53 fm for all 
nuclei, and c was determined from the known rms radii of nuclei and nucleon R rms and r 0rms 
with the help of the formula 



V' z = R z - r A = -r 2 

' ^rms ''rms rms g 



3 r^ 7/7T<i\ 2 i 
+ 3 \~c) 



(3.J 



In general, the accuracy of such a matching should be checked every time since the nuclear Fermi 
densities Pf(j) have different values of a for different nuclei. Besides, the obtained parameters 
of G-functions depend on a choice of matching points. 
Starting with [21] the symmetrized SF-function 



( ) sinh i?/a 1 1 (3 9) 

cosh R/a + cosh r/a 1 + exp 1 + exp 

has come into use, first, in calculations of nuclear form factors in eA-scattering, and next, 
in other problems of nuclear physics, too. This function has several advantages as compared 
to the F-function, allowing much room for analytical calculations j2H], !2H]- Its shape is the 
universal one for a satisfactory modelling the nuclear densities of light, middle and heavy nuclei 
[27\. It is evident from ()3.9|) that for middle and heavy nuclei (R ^> a) this function in fact 
coincides with the usual Fermi function Uf{t) = 1/(1 + exp[(r — R)/a\). Therefore, parameters 
of this function can be taken from the existing Tables of Fermi distributions for both nuclear 
densities j2Hj and point-like densities of nuclei j22]. For our task it is important that just for 
the SF-function the respective profile was obtained in an explicit form (SO] arid therefore the 
following calculations can be considerably simplified. So, the SF-density distribution and its 
profile have the following form (3.10) (3.11): 
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PsF{r) 



Psf(0) 



sinh R/a 



cosh R/a + cosh r/a 



Psf{0) 



3A 



/ira\ 2 



p SF (s) = 2R p SF (0) 



sinh R/a 



P(s) 



(3.10) 



(3-11) 



cosh R/a + cosh s/a 

Here the main dependence of the profile on s is determined by SF-function with the same 
parameters as in the density psF(r). The corrective factor P(s) is presented in [3U] and specified 
with the help of the auxiliary function x(s). This latter obeys the condition x(s) <C 1 which 
allows one to simplify P(s), so that it arrives at 



Pis) 



R 



ln(4/x(s)), 



x(s: 



cosh s/a 



K cosh s/a + cosh R/a \ cosh s/a 



1 + 



« - 1 



(3.12) 



Here k is expressed by the radius R and diffuseness a as 
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1.10315 + 0.34597(P/a) - 0.00446(P/a) 2 . (3.13) 

The numerical values of the coefficients in (J3.13)) were found in |3D1 by fitting the profile 1)3.11)1 
to numerical values of the profile integral (J2.5)) for Psf{ t ) (|3.10|) in the region of 5 < R/a < 20. 
In the center of a nucleus one has P(0) = 1, and in the region of the main contribution from 
s = R to 00 it changes a little by ~ 0A(a/R). This enables one to take P(s) at one point only, 
for example, at s = R, or at s = S1/3 = R + a In 2, where the density itself falls by three times 
6 . Then, if cosh R/a ^> k, one obtains 



p / p 

2.48945 + 0.34597 0.00446 - 

a \ a 



and, respectively, 



Psf(s) ~ 2i? p 5F (0) 



sinh i?/a 



cosh R/a + cosh s/a 



Pa(P)- 



(3.14) 



(3.15) 



In Fig.l are shown the profiles of SF-density for three nuclei: 16 0, 40 Ca, 208 P6, calcu- 
lated numerically (points) using eq. (|2.5jl . and the corresponding G-profiles (dashes), matched 
according to ()3.fij) . (|3.7j) at R and i? + 4a. Full curves represent Psf(s) calculated with the 
help of the analytic formula 1)3.11)1 . The parameters of SF-densities are taken from [2Zj and 
are given in Table 1. One can see that a strong discrepancy exists between (5-profiles and 
the initial SF-profiles in the inner region (for example, by two orders of magnitude for 208 Pb), 
and also in the region where the density falls off by two orders of magnitude and more. The 
profiles of the uniform and Gaussian forms differ much more from the profiles of the Fermi 
function. In Fig.2, the functions Psf(s) are depicted for the same nuclei but calculated within 
analytic formulae: full - by 1)3.11)) with the exact corrective factor P(s), and dashed lines - by 
eq. (j3.15j) with the approximated one P a (s = R). It turns out that the use of the corrective 
factor at the radius point practically does not change the behavior of profile functions in the 

6 It was established in ^] that the behavior of the nucleus-nucleus scattering is responsible for the region 
where densities overlap in their periphery at s > S1/3. 
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Figure 1: Thickness functions for SF-densities. Bold dots - numerical integration of eq.(2.5), 
dashes - G-profiles adjusted to the bold dots curves. Solid curves - SF-thickness densities 
calculated by analytic eq.(3.11). Parameters are in Table 1. 

peripheral region. A slight difference, no more than a factor of two for 208 Pb, appears in the 
inner region only, which is far less than it is observed in Fig.l when one uses the G(s)-functions. 



Table 1. Parameters of nuclear symmetrized Fermi density distributions psF{f, R,a) 7 . 



Nucleus 


R, fm 


a, fm 


Rrms i ftTL 


Ref. 


12 C 


2.214 


0.488 


2.496 


m 


16Q 


2.562 


0.497 


2.711 


[21\ 


20 Ne 


2.74 


0.572 


3.004 


|2B| 


24 Mg 


2.934 


0.569 


3.105 




27 Al 


3.07 


0.519 


3.06 


[28J 


28 Si 


3.085 


0.563 


3.175 




32 S 


3.255 


0.601 


3.370 




40 Ca 


3.556 


0.578 


3.493 




66 Zn 


4.340 


0.559 


3.952 


[28J 


89y 


4.86 


0.542 


4.27 


EB| 


208p b 


6.557 


0.515 


5.427 


[27] 



So, in the subsequent discussion we have a possibility of choosing the profiles p(s) of the 
explicit analytic form for different nuclear densities, namely, for the Gaussian, uniform and 
symmetrized Fermi functions. Below they will be used in calculations of microscopic eikonal 
phases x(b) and total reaction cross-sections. 

First we consider the convolution integral (|2.17|) where the thickness density of an incident 
nucleus p°(C) is taken in the Gaussian form ([3.5)1 . If one substitutes into (j2.17J) the expression 



P°G,p(\b- si) = (V^a°G, P ) Pe, P (°) ex P 



,o 2 
G,p 



(b 2 - 2bscos(f + s 2 ) 



(3.16) 



7 Parameters of the Fermi-distributions pf(t, R,a), taken from [2Hj, are very close to those of psF(r,R,a) 
for the given three nuclei 
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Figure 2: Comparison of the SF-thickness densities calculated by analytic eq.(3.11) with the 
exact correction function P(s) (solid curves), with the approximate P(s=R) by eq.(3.15). Pa- 
rameters are the same as in Fig.l. 



and integrates over tp by using the definition of the Bessel function Iq{x) [3T] of an imaginary 
argument, it follows: 

oo 

ZcA b ) = 27r (v / ^«g, p )Pg iP (0)exp(-— T J / sds p t (s) expf-—^) hf-^)- (3.17) 

G,p J Q G,p G,p 

In a more general case, when the magnitudes of rms radii of the nucleon and NN-interaction 
differ, it is convenient to use the convolution integral (J2.4)) with profiles of point-like densities 
for both nuclei. Then one can show that in the case of the G-thickness density of an incident 
nucleus f!3.5|) the convolution integral takes the form: 



a° 2 



lG,N,t{ h ) = 2 \o2 °f n 2 ^ a G,p) Ph, P (°) eX p(" 
61,=, _ tUj; \ 



I n z v ■ <^,p' ■ u^v / - \ „o 2 _i_ „2 
G,p ^ N v a G,p ^ N 



X 



/ \ / lbs \ 

sds p° t {s) expf- 2 2 ) hi 2 2 )■ (3.18) 

a G,p ^ N y a G,p~*~ a N / 

If for a target-nucleus one also takes the Gaussian function, then integration in (3.19) is per- 
formed explicitly pTT] . and one obtains [22j 



1 (v 7 ™^ ) 3 (\/7Ta%,j / /,- 

Zb^sW = ~ fl02 f a o2 + Q T P^(0)P^(0) ex P (- )■ (3.19) 

G,p + tt G,i + tt Af a G,p + tt G,t + a N / 

Note that in the case of normalized Gaussian functions (|3.1j) it is necessary to change G by G 
in (jUlD-ilH!), and put (v^Fa°) 3 p^(0) = A t . 

We imply that the realistic density distributions of middle and heavy nuclei are the (sym- 
metrized) Fermi functions. Their profile functions ()3.10j) and (J3.15|) are known in a certain 
form. The Bessel functions Io(x) and Jo(x) are also known explicitly in the form approximated 
by polynomials [S2]- So all functions in the convolution integrals ()3.17j) and ()3.18|) are given 
explicitly which is highly feasible for numerical integration. 
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Figure 3: Total reaction cross-sections for different thickness densities. The G-profile of the 
projectile nucleus 12 C is adjusted to its SF-point-like thickness density having parameters from 
Table 2. For the target-nuclei, solid curves are for the SF-densities (parameters in Table 1), 
dashes-dots - for unified densities, dashes - G-model. 

In the case when the density distributions of both the nuclei are given as the SF-functions, 
it is reasonable to use the convolution integrals in the momentum representation (|2.9|) . 1)2.14)1 
or (|2.1fij) . Then, it is convenient to take their thickness functions approximated by ()3.15j) : as 
a result, the Fourier transform can be easily calculated. Indeed, inserting ()3. 15j) into ()2.8|) one 
obtains 

p SF (k) = AtxR p SF (0) P a (R) F SF {h, a, R), (3.20) 

where 

oo 

/sinh R I a 
adaJ " ( ^ co s h fl / a + c 05Wa - < X21 > 



Taking into account the peripheral nature of nucleus-nucleus collisions one can assume that the 
main contribution is made in the region when ks ^> 1. Then (see, e.g., p3J) we have 

F SF (k) = ™ R J x {kR). (3.22) 
sinh nak 

In j2Hj, the higher order corrections to ()3.22j) are established but they do not give significant 
contributions at kR ^> 1. 

In Fig. 3 are shown the calculated total reaction cross- sect ions for collisions of 12 C with 
16 0, 40 Ca, 20S Pb at energies from 50 to 350 MeV/nucleon. For the target-nuclei, different 
nuclear densities are chosen in the form of SF-functions, modified G-functions and U-functions 
of uniform distribution. The convolution integral (j3.17|) has been used for the x _ phase. The 
profile function G for the incident nucleus 12 C was shifted at s% — c S2 = c + 4a (see (3.5)- 
(3.7)) with the profile of the corresponding point-like SF-density whose parameters are given in 
Table 2. When one tested G-profiles for the target-nuclei, the parameters of their densities psf 
were taken from Table 1. The _R rms -radii which have been used in calculations of the uniform 
distribution radius R u ()3.4|) are also quoted in Table 1. The energy-depended total nucleon- 
nucleon cross-sections a^N are taken from [23 . One can see that for a uniform density, the total 
reaction cross-sections (dot-and-dash curves) have highly different forms as compared to ones 
for the physically justified SF-densities (full lines). Besides, both calculations with the SF and 
G-models (dashed) are close to each other. A slight excess of cross-sections for SF-functions 
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Figure 4: Dependence of cross-sections on the point-like density parameters obtained from 
eA-scattering data. Solid curves - with parameters from Table 2, dashed curves - with param- 
eters given by eq.(4.4); dashes-dots - parameters from [37]; bold circles - calculations in the 
momentum representation using eq.(3.20) (see the text). 

at low energies results from their extended "tails". Instead, a weak relative increase in cross- 
sections at higher energies for the G-functions of the target-nuclei is due to their larger values 
in the inner nuclear region, which leads to the earlier "activation" of absorption as compared 
to the SF-model. 

The cross-section calculations in the momentum representation with utilizing the realistic 
/^-densities and the approximate correction P a (_R)-function (j3.14|) . and the related approxi- 
mate profiles (j3.15J) and ()3.20|) show (see Fig. 4) a slight difference ~ 2% of such calculations 
("bold" points) from the exact ones (solid curves). 

Table 2. Parameters of the symmetrized Fermi density distributions p° SF {r, c, d) 
of the point-like nucleons in nuclei [29\ . 



Nucleus 


12 C 


16Q 


24 Mg 


28 Si 


32 g 


40 Ca 


c, fm 
<i, fm 


2.275 
0.393 


2.624 
0.404 


2.984 
0.484 


3.134 
0.477 


3.291 
0.520 


3.593 
0.493 



4 Calculations and discussion 

4.1 On the use of nuclear densities from eA-scattering 

The convolution integrals ()2.4j) and ()2.9|) depend on profiles of both nuclei composed of point- 
like densities p°(s), while the transformed integrals ()2.16j) . (j2.17|) include only one density p°(s) 
whereas the other density is the nuclear one p{s). So the problem of obtaining point-like 
densities of nuclei is not elucidated completely, and we consider how they can be extracted 
from existing data. In general, one can calculate them in the framework of nuclear models. 
However, from the outset we posed the task of using primarily the data of other experiments, 
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for instance, the data on nuclear charge form factors. In this case we consider, so far, that the 
realistic point-like densities have the form p%)p( r )- 

The first method [29J is based on the representation of a nuclear ("experimental") form 
factor in the form like (|2.12|) as follows: 

F(q) = F P (q)F°(q). (4.1) 

Here F°(q) is the form factor of a nucleus with point-like nucleons and Fp(q) is the proton 
form factor presented by the dipole formulae which can be approximated at small momentum 
by the Gaussian function 

F P {q) = (l + ^% ss ) ~ 2 * exp(-g 2 r 2 _/6). (4.2) 

Then, the obtained F°(q) is analyzed within the model-independent method to obtain a point- 
like density p°(r) as a sum of the pg F (r)-function with its derivatives multiplied by the fitted 
coefficients. The latter reproduce the so-called radial variations of densities. In this procedure, 
every F°(q) with the respective trial density p°(r) is calculated in the high-energy approximation 
[34J, 35 , the analytical method which gives results in close agreement with numerical solutions 
of the Dirac equation. In Table 2 we reproduce part of data from namely, the radii c and 
diffuseness d parameters of densities p° SF (r,c,d) excluding the radial variations which play an 
important role only at large q. In j2H], the employed proton rms-radius rQ rms = 0.658 fm 2 
slightly differs from rQ rms = 0.65 fm 2 , the matter rms-radius, which was used in calculations 
of the double- folding potentials ^H]- However, by this reason the point-like densities from 
can be related to the nucleon distributions Psp{r) rather than to the proton ones 8 . 

The other method [HE] of obtaining c and d parameters of p°(r, c, d) is based on approximate 
analytic calculations of r n -moments of densities pp(r, R,a), where the latter is given by its 
folding form like (J2.ll)) . The obtained explicit results for moments are compared with those 
obtained by using the standard form Pf{t, R, a) and lead to 



R 



3V R 



d 



Y f '"O rms ^ 

~ 2 V Tra 



(4.3) 



In evaluations, terms of orders higher than (ira/R) 2 and ?"o rms /9c 2 were neglected. 
If one inserts parameters R and a of psF(r, R,a) from Table 1 and r\ rms = 0.658 fm 2 into 
equations ([4.3)1 . and compares the obtained values with the respective parameters from Table 
2, one discloses the former to be smaller as compared to those in Table 2 (about 1% for c and 
not more than 10% for d). The effect of this discrepancy on the corresponding total reaction 
cross sections is shown in Fig. 4. The solid curves are calculations with parameters c and d from 
Table 2. The dashed curves are with using c and d, estimated by eq. ()4.3|) at R and a from 
Table 1. The dash-dotted curves show the cross-sections for c and d from [37] ( 12 C: 2.1545, 
0.425; 16 0: 2.525, 0.45; 40 Ca: 3.60, 0.523), where they were used in calculations of the real 
parts of nucleus- nucleus folding potentials to explain the elastic scattering cross-sections at 
energies of about 10 MeV/nucleon. In all the cases we take r 2 Nrms = 0.658 fm 2 . One can see 
that for every set of colliding nuclei differences between the respective cross-sections occur in 

8 If one assumes the relation between proton and neutron densities to be p° N (r) — (N/Z)p° z (r) and takes 
rms-radii of densities of nuclei consisting of nucleons Rj. ms , the point-like nucleons TZ^, ms , and protons (r 2 )p = 
0.76 fm 2 and neutrons (r 2 ) N = — 0.11/m 2 (see QH), then it follows from R 2 ms = H 2 ms + (r 2 )p + (r 2 )N = 
7Z 2 ms + (r 2 ) that the rms-radius of a nucleon is (r 2 ) = (r 2 )p + (r 2 )^ = 0.65 fm 2 . 
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Figure 5: Dependence of cross-sections on the NN-interaction radius: r 2 Nrms = 0.658 fm 2 (solid 
curves), 0.630 fm 2 (dots), 0.69 fm 2 (dash-dots). Parameters c and d for 16 and Ca are in 
Table 2, and for 208 Pb they are done by eq.(4.4) with R and a from Table 1. For 12 C the 
G-profile is the same as in Fig.3. 

the limits of ~ 6 4- 10%. Nevertheless, we incline to believe that a more rigorous method of 
obtaining point-like density parameters is to analyze the F°(q) form factors of nuclei, and thus 
it is important to make up Tables of the respective densities. 

4.2 On establishing NN-interaction radius 

It was shown in Sec. 2 that at the same radii of the nucleon and NN-interaction r 2 rms = r 2 Nrms 
the convolution integral is reduced to a simpler form ()2.17j) which contains only thickness func- 
tions of nuclear matter distribution of a target-nucleus and the point-like density of a projectile 
nucleus. As to the NN-interaction parameter a 2 N = (2/3)r% rms , it is known from the scattering 
data of free nucleons to give r 2 Nrms in the limits of 0.63 4- 0.69 fm 2 . At the same time, in the 
dipole formula, the nucleon rms-radius rl rms was used as 0.658 fm 2 [22], and in calculations of 
the double-folding potential ^H] it was taken to be 0.650 fm 2 . The effect of their diffferences 
on the reaction cross-sections is seen from Fig.5. In calculations, convolution integrals were 
taken in the form ()3.18|) . The parameters c d of the point densities for 12 C, 16 0, 40 Ca are 
given in Table 2, and for 208 Pb they were evaluated using eq. ()4.3|) . For the projectile nucleus 
12 C we used G-profile adjusted to the respective thickness SF-density in the same manner as 
for Fig.3. It is seen from Fig.5 that the obtained cross-sections, in fact, coincide to each others. 
Therefore, the study of total cross-sections does not allow us to distinguish between rms radii 
of a nucleon and the NN-interaction. Moreover, one should bear in mind that the amplitude 
of scattering of free nucleons and its parameter can differ from those scattered in nuclear 
medium. 

4.3 Influence of the trajectory distortion 



In the repulsive Coulomb field the trajectory of an incident nucleus deflects from the scat- 
tering center, which results in decrease of the total reaction cross-section. This effect was taken 
into account in 38 j by replacing, in the phase x(b), the impact parameter b by the distance of 
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Figure 6: Trajectory distortion effects on cross-sections. Solid curves - no distortions; dashes 
- only the Coulomb distortion included; dashes-dots - effect of the Coulomb and nuclear dis- 
tortion. For parameters see the text. 

the turning point b c in the Coulomb field 

b -> b c = a + Va 2 + b 2 , (4.4) 

where a = Z p Z t e 2 /2E c , m , is the half-distance of the closest approach in the field Z p Z t e 2 /r 
at b = 0. Later, the procedure of exchanging b by b c was also used for the nuclear part of 
the phase $jv(6) in calculations of differential cross-sections of elastic scattering [T3], and, in 
general, it proved to be correct (see, for instance, (HH])- In addition, at peripheral collisions 
one can account for a contribution of the real part V(r) of an attractive nuclear potential, 
which brings the Coulomb trajectory closer to the target-nucleus. If one assumes the region 
b > R s = R p + Rt + (a p + a t ) In 2 to be the main for elastic scattering where the nuclear densities 
overlap less than 1/3 of their values in the center |19] . then the influence of the "tail" of the 
nuclear potential can be taken into account by exchanging 

b -> b c = a + Va 2 + b 2 , (4.5) 

where a = (Z p Z t e 2 — R s \V(R S )\) /2E c . m _. A more refined way of inclusion of nuclear distortion 
was elaborated in jlUj and applied in a series of works (see, for example, [Ej). Nevertheless, if 
the optical potential itself is obtained by numerical fitting to experimental data, then the use 
of its real part to correct the Coulomb trajectory in calculations of the reaction cross-sections 
or loses its meaning. Indeed, in these cases the data on a R are usually included into the fit 
procedure, or, if they are not available, then they themselves are calculated on the basis of the 
Sj-matrix elements obtained by fitting the differential cross-sections of elastic scattering only. 
Often these "calculated" ctr data are called "experimental" cross-sections. So the use of the 
nuclear trajectory distortion is meaningful only for construction of eikonal phases of distorted 
waves when calculating direct inelastic and nucleon removal reactions. Another situation is 
when the real part of the nucleus-nucleus potential is calculated, for example, using the double- 
folding method. Then it is reasonable to calculate both the differential and total cross-sections 
in the Glauber-Sitenko approach taking into account the trajectory distortion by both the 
Coulomb and nuclear field. 

Figure 6 demonstrates the calculations by (|3.17|) of the total cross-sections for 13 C+ 208 Pb 
and 16 0+ 28 Si without including trajectory distortion (solid curves), with the Coulomb distor- 
tion only, using the method (|4.4j) (dotted), and with the Coulomb and nuclear distortion by 
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the method (|4.5|) (dash-dotted). In the last case one should specify the parameters of the 
nucleus-nucleus potentials, and for 13 C+ 208 Pb we took them from ^1] (the " C-potential" at 
390 MeV), and for 16 0+ 28 Si from 42 j (the " E-potential" at 215.2 MeV). In computing phases, 
we used for 13 C the parameters of the 12 C point-like density from Table 2. The parameters 
of 16 are also given in Table 2. For 208 P5 and 28 Si, the parameters were computed with the 
help of (|4.3|) using the values R and a from Table 1. As expected, the inclusion of Coulomb 
distortion leads to appreciable corrections, of the order of 10100 MeV/nucleon and less, which 
is beyond experimental errors. As the collision energy increases, these corrections diminish and 
for lighter targets- nuclei (the right-hand part of Fig. 6) they reach about 2they can be neglected. 
Contribution of the nuclear distortion is poor in comparison with the Coulomb one for 208 Pfr 
and comparable with this for the reaction on the 28 Sz-target, although in the latter case both 
the effects give small contributions. Besides, it is necessary to bear in mind that the real part 
of the nuclear potential decreases with increasing energy, but this was not taken into account. 

4.4 In-medium effects 

In the microscopic approach we deal with the total cross-section cr NN (ei a b) of the NN- 
scattering of free nucleons. This cross-section depends on energy, and thus defines the main 
dependence of the nucleus-nucleus cross-section on the collision energy Eiab = ei a bA p . We have 
taken the parametrization of (JNN^iab) from j2H| in the energy interval = 10 MeV 4- 1 GeV . 
More generally, one should take into account the in-medium effect on the nucleon-nucleon 
interactions in nuclear matter. Usually, for this aim the cross-section is multiplied by the 
factor f m , and therefore, in the convolution integral the isospin averaged cross-section ()2.3|) is 
exchanged as follows: 

<?np => CF np X f m (np), x f m {nn). (4.6) 

The factors f m (np) and f m (nn) depend of the nucleon energy ei a b = Eiab/A and on the density 
of nuclear matter. The problem of the in-medium corrections of the NN-forces has been inves- 
tigated in many works. So in |43 s , based on the Dirac-Bruckner approach for nuclear matter, 
numerical calculations were made of the total NN cross-sections. Parameterization of these 
calculations was given in jH] in analytic form for the correction factors 

1 _l_ 9fl 88 c- 04 n 2 - 02 1 _L 7 779 c-0.06 -1.48 

f , ) = 1 + 20.88 p 1 + 7.772 e lab p 

Jm[n P ) i + 35. 86 p i.9o > M nn ) 1 + 18.01 p 1 - 46 ' 1 ' 

Here the energy is given in MeV, and the density in fm~ 3 . One can see that for free nucleons 
when p = 0, we have f m ( n p) = fm( nn ) — 1> an d as the density increases these factors f m vanish, 
as well as the corresponding effective cross-sections 9 . It is difficult to compute the convolution 
integrals with correction factors in the form (J4.7)) including the dependence on r, and we limit 
ourselves to qualitative estimations only. In Fig. 7 are shown such calculations when the values 
of nuclear densities in ()4.7|) are assumed to be constants p = p = p p + p t for every region where 
the colliding nucleons can be. Then, denoting by p = p p (0) + Pt(0) the net density in the 
centers of colliding nuclei, we show the calculated cross-sections when the f m factor are taken 

9 In the microscopic models of the double-folding potentials the same problems arise when constructing the 
energy and in-medium dependence of NN-potentials. However, the advantage of the Glauber-Sitenko approach 
is that the main dependence of energy is already included in parametrization of the total cfnn cross-sections of 
free nucleons. 
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Figure 7: In- medium effect on cross-sections. Solid curves - no in-medium dependence (p=0); 
bold dots are when the density in the overlap region is taken as l/20p o with p = p p (0) + pt(0) 
expressed through densities at the centers of nuclei; dots - for overlap with l/6p ; dashes-dots 
- for l/3p ; dashes - at p . 

for free nucleons (at p=0, solid lines), and also for nucleons in the medium at p=(l/20)p o (bold 
lines), (l/6)p Q ) (dots), (l/3)p D (dash-dotted lines), p Q (dashed curves). It is seen that the 
inclusion of the in-medium factors f m can diminish the total reaction cross-section by 4-7% and 
that the dependence on the density turns out to be strongly nonlinear. 

5 Conclusions, comparison with experimental data 

1. It is shown that the initial expression for the microscopic eikonal phase can be represented 
in a convenient form where one of two density thickness functions of point-like nucleons in nuclei 
transforms into the profile of a nuclear density obtained from independent experimental data, 
e.g., from nuclear form factors. Ultimately, it turns out that in calculations of nucleus-nucleus 
cross-sections it is not necessary to introduce free parameters but, instead, to base on the known 
experimental data, only. 

2. The main problem of the microscopic approach is to calculate eikonal scattering phases. 
For this aim, many people use the simplest nuclear densities in the form of Gaussian functions. 
Instead, here we demonstrate how one can use the realistic density distributions in the form 
of (symmetrized) Fermi functions, whose parameters are known for many nuclei. In this case, 
all functions in the eikonal integrals turn out to be given in the explicit form, which highly 
simplifies numerical calculations. 

3. Two methods of obtaining the radius c and diffuseness d parameters of the point-like 
density p SF {j^ c, d) are presented. One of them deals with nuclear form factors, and the other 
obtains these parameters from the known nuclear densities psf{t,R,cl)- It turns out that 
the difference between the cross-sections calculated by using these two sets of parameters can 
achieve the value beyond the bars of typical experimental cross-sections. We conclude that 
a more justified method is that based on the analysis of form factors of nuclei obtained by 
subtracting the nucleon form factor and, if necessary, the form factor of the nuclear center-of- 
mass motion, as was carried out, for example, in [29]. 

4. In calculation, it was established that the rms radii of the nucleon-nucleon interaction 
r Nrms an d °f the nucleon itself rQ rms may be considered to be equal. The slight difference 
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between them is in the range of accuracy of their evaluation, and this does not affect practically 
the results of calculations of the nucleus-nucleus total reaction cross-sections. At the same time 
& t r Nrms = r orms the convolution integrals take highly simple forms (|2.16|) . (|2.17|) where only 
two profile functions overlap, one is for the nuclear density and the other is for the point-like 
density. Next, the convolution integral takes a simple one-dimensional form. 

5. In many typical cases one should take into account a Coulomb trajectory distortion by 
means of the formal replacing, in the phase x{b), the impact parameter b by b c , according to 
(|4.4jl . Additional distortion of the trajectory by the tail of a nuclear potential is not reasonable 
because its parameters are to be fitted, in particular, to the same data on total cross-sections. 

6. The question whether the in-medium factor f m for corrections of the NN-cross-sections 
should be taken into consideration remains open for us till now. The estimations show that 
at intermediate energies this factor does not change substantially the nucleus-nucleus total 
cross-sections. Nevertheless, the use of the above given factors f m (np) and f m (nn) from |H] 
hinder the usual calculations. Other authors faced similar obstacles, for example, when solving 
a simpler problem of the pA scattering data analysis jlH]- Besides, in both these works, the 
Gaussian functions were chosen as the basic nuclear densities, while their replacing by the 
realistic ones can change the conclusions about the structure of the factors f m as compared to 
those given by eq. ([4.7)1 . For example, successful agreement of calculations with the data on 
total reaction cross-sections for a series of reactions pA, a A and 12 C + 12 C was obtained in 
with realistic densities and with the factors f m differing in form from those given above 10 . 

7. Our calculations show that in the case of relatively light incident nuclei, it is more 
profitable to use the convolution integral ()3.17[) with the thickness function p°(s) for the point- 
like density of the projectile nucleus in the form of the modified Gaussian G-function. The latter 
is determined with the help of ([3.6)) and ()3.7|) by means of parameters of realistic SF-function 
taken, for example, from Table 2. In the case of heavier incident nuclei it is advantageous 
to use the convolution integral ([2.16)1 where one can insert the known explicit form of the 
Fourier-Bessel profiles of realistic SF-densities (|3.22|) for both the nuclei. In all calculations the 
Coulomb distortion of trajectory should be first of all taken into account. 

8. We emphasize that in the existing literature no attention is focused on the problem of 
correct use of the initial formula for the convolution integral. From the above discussion it is 
evident that if the finite radius of the NN-interaction (a 2 N ^ 0) is taken into account explicitly, 
both the densities must be taken as densities of nuclei for the point-like nucleons. Further, 
when the rms-radii of the nucleon and the NN-interaction are equal to each other, the NN- 
interaction factor disappears in the convolution integral, but one of the point-like densities is 
transformed to the nuclear one. However, in this case, the absence of the NN-factor /(£) does 
not mean that one uses the zero-range approximation. The confusion arises also, when one 
assumes a 2 N = (2/3)r^ rms = in ([2.4)1 and calls this case the zero-range approximation, while 
at the same time takes for both the densities the table data, i.e., just the nuclear densities, 
instead of the point-like ones. 

In Fig.8a,b, our calculations are compared with the experimental data collected in [49]. The 
SF-density parameters are given in Table 1. For the density p° SF the parameters c and d of 
the projectile 12 C are given in Table 2, and for the incident nucleus 20 Ne they are calculated 
with the help of ([4.3)1 and the data from Table 1. The convolution integral has been taken 
in the form of ([3.17)1 . Only the Coulomb trajectory distortion was taken into account. Thus, 

10 In calculations of the real part of the double-folding nucleus-nucleus potential, the in-medium effect on 
the NN-potential is parametrized by simpler factors f m in the form of the step and exponential functions of 
densities p(r) (see., e.g., 07]). 
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Figure 8: Calculated cross-sections compared with experimental data from [48]. Parameters 
are from the eA-scattering data (see in the text). Only the Coulomb distortion is included. 
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no free parameters have been introduced. One can see that in all cases a good agreement 
with experimental data is observed, except for 12 C + 89 Y and 20 Ne + 27 Al. In the latter 
cases, some discrepancies can be due to the fact that the parameters for densities of the odd 
nuclei 89 Y and 27 Al were obtained from eA form factors (for references see [28]) with the help of 
formulae for spinless nuclei. However, if one introduces the in- medium factor f m , the mentioned 
disagreements can be removed. Nevertheless, we believe that first of all one should improve 
the data on geometric parameters of the given nuclei. The other remark is that calculations 
of cross-sections in |3H] with the help of Gaussian functions which reproduce the behavior 
of "tails" of densities, give enhanced values which are beyond the experimental bars shown 
in Figs.8a,b. This is conceivable that in [1H] they took the nuclear densities instead of the 
point-like ones, as it was mentioned above. On the other hand, calculations in [H| with the 
uniform density distributions for the target-nuclei and Gaussian functions for the projectiles 
give underestimations of the cross-sections for 20 Ne + l2 C and 12 C + 27 Al, in spite of that they, 
probably, use not the point-like but nuclear densities for estimations of rms-radii needed to get 
the step radius R u in (|3.3j) . So, this result confirms that shown in Fig. 3. 

In conclusion, the authors would like to thank the Infeld-Bogolubov Foundation for its 
support, and E.B.Z. acknowledges the support of the Russian Foundation for Basic Research 
(grant 0001-006-17). 
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